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ABSTRACT 

Magnetic fields are important at every scale in the star formation process: from the dynamics of the ISM in 
galaxies, to the collapse of turbulent molecular clouds to form stars and in the fragmentation of individual star 
forming cores. The recent development of a robust algorithm for MHD in the Smoothed Particle Hydrodynamics 
method has enabled us to perform simulations of star formation including magnetic fields at each of these scales. 
This paper focusscs on three questions in particular: What is the effect of magnetic fields on fragmentation 
in star forming cores? How do magnetic fields affect the collapse of turbulent molecular clouds to form stars? 
and: What effect do magnetic fields have on the dynamics of the interstellar medium? 

Key Words: Stars: formation — Stars: Pre-main sequence — magnetohydrodynamics (MHD) — galaxies: magnetic 
fields — galaxies: spiral 



1. INTRODUCTION 

Magnetic fields are observed to be important at 
every scale in the star formation process-in the dy- 
namics of the ISM in spiral galaxies, during the col- 
lapse of individual molecular clouds to form stars, 
and even down to the fragmentation scale in star 
forming cores. However we remain significantly ig- 
norant of the role of magnetic fields in constraining 
the dynamics of collapse at each of these scales. 

In this contribution I will discuss our efforts to 
incorporate magnetic fields into simulations of the 
star formation process at these three different scales: 
in star forming cores, in collapsing turbulent molec- 
ular clouds, and in the dynamics of galactic spiral 
arms. These can be phrased into three unanswered 
questions regarding the role of magnetic fields in star 
formation (which we will discuss in order of smalls 
large scales): 

1 . What is the effect of magnetic fields on fragmen- 
tation in star forming cores? 

2. How do magnetic fields affect the collapse of 
molecular clouds to form stars? 

3. What effect do magnetic fields have on the dy- 
namics of the interstellar medium? 

We have been able to start addressing these three 
questions in our simulations because of the develop- 
ment over the last 6 years or so (from my [Price's] 



PhD thesis onwards) of a workable algorithm for in- 
corporating magnetic fields into the Smoothed Parti- 
cle Hydrod ynamics (SPH) method. SPH (for recent 
reviews see IPricei l2004t iMonaghanl 120051 ) is ideally 
suited to simulations of star formation because of its 
Lagrangian nature which means that resolution nat- 
urally follows mass. Thus, given that the primary 
computational challenge of star formation is deal- 
ing with the huge range in length and time scales 
involved during the collapse of a diffuse molecular 
cloud to form a star, SPH makes a natural choice, 
though attempts to add magnetic fields have had a 
long and somewhat tortuous history. 

2. SMOOTHED PARTICLE 
MAGNETOHYDRODYNAMICS 

In principle an algorithm for Smoothed Par- 
ticle Magnetohydrodynamics can be derived 
straightforwardly from a variational principle 
( Price fc Monaghanll2004rJ ) star ting with the usual 
hydro dynamic Lagrangian ( Monagh an fc Pried 
1 2 lh with the magnetic energy subtracted: 
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Constraining the Lagrangian using the SPH density 
summation, the first law of thermodynamics and the 
SPH representation of the induct ion equation for the 
magne tic field leads naturally (jPrice fc Monaghan 
2004bJ) to the equations of motion for the magnetic 
fields in SPH in a form which conserves momentum 
and energy exactly (in practise this means to the 
accuracy of the time-stepping algorithm). This is 
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of course very beautiful and elegant since it is rare 
to be able to derive anything, let alone the equa- 
tions of motion themselves, from a first-principles ap- 
proach in any numerical method. Unfortunately for 
life (mine especially) and probably related to some 
deep conservation law regarding the conservation of 
beauty (perhaps the divergence-free nature of ele- 
gance?), that is where the beauty ends and the fun 
begins, as there are a number of thorny technical is- 
sues which must be dealt with from that point before 
the algorithm can be deemed 'workable'. 

Briefly, these issues are that: i) the exactly- 
momentum conserving representation of the Lorentz 
force in SPH turns out to be numerically unstable 
in the regime where magnetic pressure dominates 
over gas pressure (the solution is to use an alter- 
native representation which is stable in all regimes 
but compromises momentum conservation slightly); 
ii) shocks in MHD are much more complicated than 
their hydrodynamic counterparts and require careful 
treatment (the solution is to apply artificial dissipa- 
tion terms including an artificial viscosity, resistiv- 
ity and possibly conductivity where appropriate) ; iii) 
that the smoothing length in SPH and MHD-SPH is 
a spatially variable quantity and thus terms involv- 
ing the derivative of the smoothing length should 
be accounted for in order to retain the Hamilto- 
nian nature of SPH (ie. exact conservation of en- 
ergy and entropy). The most thorny of all the tech- 
nical issues (perhaps the rose bush itself) is that 
nasty fourth Maxwell equation V • B — which 
blights all numerical MHD by the simple fact that it 
enters the MHD equations as a constraint on the 
initial condition which should remain satisfied for 
all time. Whilst divergence cleaning methods sim- 
ilar to those adopted by several grid-based MHD 
schemes can also be app lied in an SPH context (see 
Price fc Monaghanll2005l for details), we have found 
a more effective approach is to use prevention rather 
than cure by formulating the magnetic field such 
that it is divergence-free by construction. A natu- 
ral choice for a Lagrangian code is to use the 'Euler 
potentials' 

B = Vax V/3, (2) 

which is divergence free by construction and has the 
advantage that the induction equation expressed in 
terms of the Euler potentials takes the form 



da d(3 

~dt ~ '~dt 



0; 



(3) 



to a vector potential formulation of the magnetic 
field where a = A z and (3 = z. 

Whilst the Euler potentials formulation is a nat- 
ural and elegant choice, the disadvantage is that in 
our simulations field growth is lost once there is no 
longer a clear mapping from the particle distribution 
at t = to the present simulation time. For ex- 
ample, imagining two rings of particles carrying the 
Euler potential representation of a magnetic field, it 
is obvious that rotating the inner ring through one 
complete rotation will return the magnetic field to its 
initial configuration (because the fields have simply 
been advected), whereas clearly in reality the field 
should 'remember' the winding. Thus field growth 
can be underestimated when using the Euler poten- 
tials formulation. 

The SPMHD algorithm using both the Eu- 
ler potentials and divergence cleaning methods 
have been extensively tested against a wide 
range of problems used to ben chmark many grid- 
based MHD codes in one (Price fc Monaghan 
2004alfbl: IPricd I2004J: IRosswog fc Pried 120071) . two 



(|Price fc Monaghanl l2005t IRosswog fc Pried 120071 ) 
and three (JRosswog fc Pried 120071 ) spatial dimen- 



corresponding to the advection of magnetic field lines 
by Lagrangian particles. Note that in two dimen- 
sions the Euler potentials formulation is equivalent 



3. THE EFFECT OF MAGNETIC FIELDS ON 
FRAGMENTATION 

One of the first questions we have looked at is: 
What is the effect of magnetic fields on fragmen- 
tation? Whilst some authors have suggested that 
magnetic fields tend to supp ress fragmentation (e.g. 
Hosking fc Whitworthl 12004 ) others have suggested 
that the presence of magn etic fields m ay in fact en- 
hance fragmentation (e.g. lBosgl2002n . To this end 
we have performed simulations of single and binary 
star formation in t he presence of a magnetic field 
( Price fc Batdl2007l ) starting from a uniform density 
1MO dense core (R = 4x 10 16 cm)(2674 AU) embed- 
ded in a warm, low density medium. The core itself 
has an initial temperature of 10K, is in solid body 
rotation and is threaded by an initially uniform mag- 
netic field either perpendicular to or parallel to the 
rotation axis. We employ a barytropic equation of 
state which is isothermal until the density reaches 
10~ 14 g cm" 3 where it becomes polytropic (7 = 7/5) 
in order to mimic the effects of radiative transfer. 
A binary star formation simulation is initiated by 
applying a 10% m = 2 density perturbation to the 
otherwise uniform density core. 

The results of a series of calculations with ini- 
tial magnetic field parallel to the plane of rotation 
(a B x field in our simulation coordinates) are shown 
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Fig. 1. Effect of magnetic fields on fragmentation: simulations of binary star formation using an initial magnetic field 
perpendicular to the rotation axis. 



in Figure [T] Wc find that in general fragmentation is 
dramatically suppressed in the presence of even rel- 
atively weak magnetic fields (mass-to-flux ratios of 
10 or lower), leading to single star formation where 
a binary would have formed in the hydrodynamic 
case. In these calculations where the initial core ro- 
tation is relatively low, the effect is mostly driven 
by the increase in the effective thermal energy of the 
cloud (ie. the extra pressure supplied by the mag- 
netic field) . This is evident for example if the calcu- 
lations arc re-run using only the magnetic pressure 
component of the force where a similar suppression 
of fragmentation is found. Actually the picture is 
slightly more complicated than that, as may be ob- 
served in Figure [1] At t/tft = 1.36 (fourth column 
in Figure [1]) the binary separation may be observed 
to decrease with increasing magnetic field strength 
until no binary is formed at all. In simulations using 
a field parallel to the rotation axis (ie. a B z field) 
this suppression of fragmentation (ie. the formation 
of only a single star) occurs for even weaker field 



strengths. However with the field in the plane of ro- 
tation as in Figure [1] the binary is prevented from 
merging by a "magnetic cushion" which is formed 
between the two overdense regions by the wrapping 
up of the magnetic field lines. Thus in this case the 
tension force actually helps to promote fragmenta- 
tion by preventing overdense regions from mergin g 
- a picture advocated previously bv IBossI (|2002l ). 
The difference is that although we indeed find that 
tension can dilute the suppression of fragmentation 
due to magnetic pressure, the net effect is always 
that fragmentation is suppressed by the magnetic 
field relative to hydrodynamics. 

We find a similarly dramatic effect on the forma- 
tion of circumstcllar discs in the presence of mag- 
netic fields - in fact disc formation can be almost 
completely suppressed even at relatively low mag- 
net ic field strengths ( mass-to-flux ratios of < 5 - 
see Price fc Bate! 120071 for more details). IShu et al. 



(2006) have suggested that this may simply be 
a manifestation of the well-known "magnetic flux 
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problem" in star formation and therefore that to pro- 
duce discs one must invoke non-ideal MHD processes 
(specifically Ohmic resistivity) to dissipate the flux. 

4. THE EFFECT OF MAGNETIC FIELDS ON 
STAR CLUSTER FORMATION 

Given that magnetic fields can have somewhat 
dramatic effects on the fragmentation of cores, what 
are the effects on larger scales'' In particular how 
does the presence of a magnetic field affect the de- 
cay of turbulence in molecular clouds? How do mag- 
netic fields affect the star formation rate and/or effi- 
ciency? What is the effect on the Initial Mass Func- 
tion from the presence of a magnetic field on larger 
scales? From a theoretical point of view, how do 
magnetic fields change the extre mely dynamic pic - 
ture of star formation painted by iBate et al.l (J2003I ) , 
where stars compete with each other for accretion 
from the parent cloud? 

4.1. Key parameters 

The importance of magnetic fields to this prob- 
lem can be quantified in terms of three key dimen- 
sionless parameters. These are: 

1. The mass-to-magnetic flux ratio, expressed in 
units of the critical value 



M 



(H\ 
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(4) 



which expresses the competition between mag- 
netic fields and gravity (sub-critical mass-to-flux 
ratios are stable against collapse, whereas super- 
critical ratios are unstable). 

2. The plasma (3 (ratio of gas-to-magnetic pres- 
sure), 







C 2 S P 



B72/V 



(5) 



which expresses the relative importance of mag- 
netic fields as a source of pressure. 

3. The ratio of RMS turbulent velocities to the 
Alfvcn speed 



/, 



turb 



0~RMS 
VA 



(G) 



expressing the importance of magnetic fields 
with respect to turbulence (magnetic fields con- 
trol turbulence when sub-Alfvenic, whereas tur- 
bulence controls magnetic fields when super- 
Alfvenic) . 



What is not generally appreciated is that these 
three parameters can be determined independently 
of each other, since, though they all have a depen- 
dance on the magnetic field strength, they depend on 
this strength to different powers (mass-to-flux ratio 
depends linearly on field strength, whereas the other 
parameters depend on the field strength squared) 
and also on different hydrodynamic parameters (for 
example changing the sound speed changes only (3 
whilst leaving the other parameters unchanged, simi- 
larly changing the turbulent RMS velocity amplitude 
changes only fturb)- The implication is that mag- 
netic fields can be dominant in one or more senses 
but not necessarily in the others. For example a 
cloud can be supercritical (magnetic fields not pre- 
venting collapse) and super- Alfvenic (magnetic fields 
not dominating turbulence) but be in the regime 
where (3 < 1 and thus have magnetic fields acting as 
the dominant source of pressure. Indeed this regime 
is perhaps the most interesting as observationally it 
i s the one in which molecular clouds a ppear to lie 
(|Crutchedll999t iHeiles fc Trolandll2004[ ). 



4.2. Initial conditions 

In order to study the effect of magnetic fields 
on star cluster formation we have performed sim- 
ulations starting f rom the same hyd rodynamic ini- 



as 



Bate et all (J2003h -- that 



is an 



tial conditions 
isolated 5OM spherical cloud of diameter 0.375pc 
(giving nn 2 = 3.7 x I0 4 ) with a temperature of 
I0K and an imposed turbulent velocity field with 
power spectrum P(k) ex fc -4 and RMS Mach num- 
ber 6.7. The equatio n of state is again barytropic 
as in lBate et al.ll2003l ). which is isothermal down to 
a density of 10 _13 g cm -3 and then becomes poly- 
tropic with 7 then changing with density in order to 
parameterise the effects of radiative transfer at high 
densities where the gas becomes optically thick to 
radiation. 

The magnetic field is imposed as an initially uni- 
form field, with strength parameterised in terms of 
the mass-to-flux ratio in units of the critical value. 
We have performed a series of 5 calculations us- 
ing M/$ = oo, 20, 10, 5 and 3 (the calculations are 
performed in ideal MHD so runs with subcritical 
mass-to-flux ratios would not collapse — rather they 
bounce and then expand freely). In all cases the tur- 
bulent velocities are super- Alvcnic (by factors rang- 
ing from oo down to 2 for the 5 runs) meaning that 
magnetic fields do not control the turbulent veloc- 
ity field. However of the 5 runs, the two strongest 
field runs (7\//$ = 5 and M/$ = 3) have initial 
plasma /3's of 0.7 and 0.25 respectively, ie. they lie 
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Fig. 2. Effect of magnetic fieids on star ciuster formation: column density in the cloud is shown for the sequence of 5 
calculations of progressively increasing magnetic field strength (left to right) at intervals of 0.2 free-fall times (top to 
bottom). Note in particular the "stripiness" in the lower parts of the cloud in the strongest field runs at t ~ 0.5 free-fall 
times and the large-scale voids evident at t/tff > 1. 



in the regime where magnetic pressure is dominant 
over gas pressure. 



4.3. Results 

4.3.1. Cloud structure 

Results from the 5 simulations are shown in Fig- 
ure [U showing snapshots of column density at in- 
tervals of 0.2 free-fall times (top to bottom) for the 
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runs of progressively increasing field strength (left to 
right). The initial magnetic field is oriented paral- 
lel to the y-axis in this figure. At first glance the 
global structure of the cloud is similar in all of the 
runs, which is a result of the fact that gravity rather 
than magnetic fields plays the dominant role in de- 
termining the initial cloud structure. Some small 
differences in the development of the initial density 
fluctuations are visible at t/tff = 0.2 - namely that 
the stronger magnetic field runs tend to produce 
smoother structures due to the additional pressure 
in the cloud. At later times, however, there are two 
particularly striking effects which only occur in the 
two strongest field runs (ie. those with j3 < 1). 

4.3.2. Magnetically supported voids 

Perhaps the most striking features are the large 
voids visible in the cloud structure at t > 0.8i// in 
the strong (3 < 1 magnetic field runs which arc com- 
pletely absent from the hydrodynamic and weak field 
runs. These voids arc directly related to the addi- 
tional support provided by the magnetic field to the 
low density material in the cloud, shown in more de- 
tail in Figure [3] which shows a close up of the void 
structure in the M/$ = 5 run at t/tff = 1.05. 
The plot shows the column density in the cloud 
(top) and the integrated magnetic pressure (bot- 
tom) which appears almost as the inverse of the 
top panel. The physics behind the creation of such 
"magnetically-supported voids" in the simulation is 
relatively straightforward: since gas is tied to the 
magnetic field lines, gas motions (due to the pres- 
ence of a turbulent velocity field) are channelled 
by the magnetic field to form dense filaments (with 
relatively higher mass-to-flux ratios and higher /3), 
leaving behind regions which are relatively lower in 
density but retain the same field strength. Thus 
in these regions the gas pressure can be strongly 
reduced whilst the magnetic pressure remains un- 
changed, producing an enhancement relative to the 
gas pressure exactly as observed in Figure [3l The 
creation of the void structure is perhaps better ap- 
preciated by viewing an animation of the magnetic 
field structure evolution overlaid on the column den- 
sity plots which is a available from the authors' web 
pag<E 

Thus the appearance of magnetically supported 
voids in the cloud structure is an inevitable and 
generic consequence of turbulence in the presence 
of magnetic fields which are strong enough to resist 
compression by the gas motion. 




2 http://www.astro.ex. ac.uk/pcoplo/dpricc/pubs/mclustcr/ 



Fig. 3. Close-up of the magnetically supported void in 
the M/$ = 5 run, showing column density (top) and 
integrated magnetic pressure (bottom). 

4.3.3. Column density striations 

The second effect which the magnetic field im- 
prints on the cloud structure in the (3 < 1 regime is 
visible in the M/$ = 5 and 71//$ = 3runs at around 
t ~ 0.5 free-fall times as a "stripy" structure in the 
low density outer regions of the cloud (visible in the 
lower part of the t/tff = 0.4 and t/tff = 0.6 panels 
in Figurc[5|). These "column density striations" are a 
result of streaming motions of the gas along the mag- 
netic field lines. This is nothing more than the well- 
known anisotropy of turbulent motions in the pres- 
ence o f a strong magnetic fi eld iGoldreich fc Sridhar 
(|l995[ k lVestuto et al.l (J2003l ) which only becomes vis- 
ible once the magnetic fields are strong enough to 
resist distortion by the gas flow and thus channel 
gaseous material in the direction of the ordered large 
scale field. 

4.3.4. Effect on the Initial Mass Function 

The additional support given to the cloud by the 
magnetic field also has a dramatic effect on the star 
formation rate, even for relatively weak magnetic 
fields. For example by t = 1.5 free-fall times around 
the total mass in stars in the M/$ = 5 and M/$ = 3 
runs is ~ 4M Q and ~ 2M Q respectively compared 
to ~ 8Mq in the hydrodynamic run. Though we do 
not follow the evolution of the magnetic field deep 
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TABLE 1 

RATIO OF BROWN DWARFS TO STARS 
FORMED IN EACH OF THE SIMULATIONS. 



into the cores at the resolution of the global simu- 
lations due to a somewhat high numerical resistiv- 
ity (and thus may underestimate the effect of the 
magnetic field on fragmentation), it is nonetheless 
interesting to examine what effect the magnetic field 
support given to the cloud on larger scales has on 
the mass distribution of stars which form. Given 
that the cloud only contains 50M Q of material in 
the first place, it is not possible to provide meaning- 
ful statistics on the detailed form of the initial mass 
function - though an effect can still be discerned by 
looking at cruder statistics such as the ratio of stars 
to brown dwarfs. 

To this end the number of stars and brown dwarfs 
and the ratio between the two is given in Table [T] for 
the 5 calculations. We find that in general the pres- 
ence of the magnetic field tends to reduce the number 
of low mass objects formed in the simulations. We 
attribute this effect to the more quiescent nature of 
the star formation in the stronger field runs, lead- 
ing to fewer multiple systems and thus fewer ejec- 
tions (which in the competitive accretion model cre- 
ate the brown dwarfs by preventing any further ac- 
cretion onto them). This in much better ag reement 
with observations fo r the substellar IMF ( Kroupal 
2001[ IChabrierll2003h compared to the previous (and 
more recent high resolution) hydrodynamic calcula- 
tions which are now known to produce a statistically 
significant excess of brown dwarfs with respect to 
observations (Bate, 2008, private communication). 



4.3.5. Comparisons with Taurus 

We have found that, even in magnetically su- 
percritical clouds, magnetic fields can still provide 
a dominant role in determining cloud structure in 
the regime where 8 < 1. The question to ask is 
then: What parameter regime is most realistic for 
the molecular clouds in our own Galaxy in which star 
formation is observed to occur? Though the deter- 
mination of sub-critical vs supercritical and super- 
Alvenic vs sub-Alfvenic is somewhat controversial on 



the cloud scale, measurements of magnetic fields in 
molecular clouds almost without exceptio n indicate 
that they lie in the regime where 8 < 1 (ICrutcher 
19991 : iBourke et al.ll200lHHeiles fc Trolandll2004h . If 



this is indeed the case we should therefore expect to 
find exactly the kind of effects on the cloud structure 
due to a magnetic field which dominates over the gas 
pressure that we find in the simulations. 

It is with this in mind that we turn to rc- 
cent observation s of th e Taurus molecular cloud by 
~ ( 20081) . 



Goldsmith et al 



Two features highlighted 
by the authors of the survey stand out: striations in 
the gas column density observed in the 13 CO maps 
which are extremely well correlated with the direc- 
tion of the large scale ordered magnetic field and 
large scale voids (or "cavities") in the cloud struc- 
ture. The presence of both of the features we find 
in the simulations in the 8 < 1 regime strongly 
suggests that star formation in Taurus is taking 
place in a regime where magnetic pressures are dom- 
inant over gas pressure. This is in broad agreement 
with direct Zeeman-splitting measu rements of field 
strengths from dense cores in Taurus (|Crutcherlll999l ) 
giving lower limits of 8 > 0.06. 

With regards to magnetically-supported voids, 
in an earlier conference proc eedings the same au- 



thors Goldsmith et al. 



( 20051 ) comment on a "very 
interesting feature" at 4 /l 30 m + 25 dcg which "ap- 
pears as a hole" , where it appears that "some agent 
has been responsible for dispersing the molecular 
gas". We suggest that this is a region where mag- 
netic pressure is enhanced relative to gas pressure 
and that the field geometry should be reasonably 
ordered with a strong line-of-sight component. The 
only maps which give any indication of such a feature 
arc the somewhat controversial observations of "fara- 



day scr eens" towards Taurus b y IWolleben fc Reich 
(l2004ah (see IWolleben fc Reichl (|2004bJ) for the de- 
tailed Taurus maps) which indeed shown a tantalis- 
ing enhancement in emission in this area. However 
more observations are needed to assess our predic- 
tion. 

5. THE EFFECT OF MAGNETIC FIELDS ON 
THE DYNAMICS OF THE INTERSTELLAR 
MEDIUM 

If we arc to ever fully understand star formation 
in molecular clouds then we must make the link to 
the larger scales in terms of the dynamics of the spi- 
ral arms behind which the molecular clouds them- 
selves are formed. On galactic scales magnetic fields 
also play a significant role in the dynamics of the 
interstellar medium and it is with this in mind that 
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Fig. 4. Effect of magnetic fieids on the ISM: gas and magnetic field structure in a two-phase ISM simulation using 
relatively weak (left) and strong (right) magnetic fields. Magnetic fields suppress structure produced by the cold phase 
of the ISM but do not eliminate it, whilst at weaker strengths the gas is able to induce significantly more disorder in 
the magnetic field. 

we have been incorporating magnetic fields into sim- 
ulations of the ISM on these scales in order to as- 
sess both the effect of the magnetic field on the gas 
flow and vice- versa (jDobbs fc PriceH2008l ) . results of 
which arc shown in Figure |U We find that the inclu- 
sion of a cold phase is crucial in performing simula- 
tions of the ISM in galaxies as even relatively small 
amounts of cold gas produce significant velocity dis- 
persion in the gas and correspondingly induce dis- 
order in an otherwise well-ordered global magnetic 
field. In general the presence of a magnetic field 
tends to inhibit the formation of structure in the disc 
(Figured]), though provided a cold phase is included 
in the calculations spurs and spiral arms clumps re- 
main present provided /3 > 0.1 in the cold gas. 

The link to observations can be made by produc- 
ing synthetic polarisation maps based on the mag- 
netic field and gas distributions in our simulations. 
An example is shown in Figure [5] which shows con- 
tours of the total synchrotron intensity and B- vectors 
for a representative case. Whilst the enhancement in 
total synchrotron emission in the spiral arm agrees 
with observations of magnetic fields in spiral galax- 
ies, we also find a corresponding enhancement in the 
ordered component of the field (though production 
of ordered fields in the interarm regions is lost when 
using the Euler potentials formulation) which is not 
observed, suggesting that further work is needed to 
understand the interaction between spiral arms and 
magnetic fields in the ISM. 




Fig. 5. Synthetic polarization map from one of our sim- 
ulations showing total synchrotron intensity (contours) 
and the polarisation B-vectors (lines) superposed on a 
column density map. 



6. SUMMARY 

What is the effect of magnetic fields on frag- 
mentation? We have found that the net effect of 



the magnetic field is always to suppress fragmenta- 
tion, driven by the extra pressure supplied by the 
magnetic field in a collapsing core. However, mag- 
netic tension can act to dilute this effect to some ex- 
tent, depending on the field geometry, by preventing 
over dense regions from merging (an effect we refer 
to as "magnetic cushioning"). Magnetic fields also 
have a dramatic effect on circumstellar disc forma- 
tion, though we caution that this may be an arti- 
fact of employing ideal-MHD calculations in a regime 
where non-ideal effects are almost certainly impor- 
tant. 
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How do magnetic fields affect the collapse of 
molecular clouds to form stars? Even mag- 
netic fields which do not prevent clouds from collaps- 
ing can have a significant effect on the cloud struc- 
ture in the regime where magnetic pressure domi- 
nates over gas pressure (where molecular clouds are 
ubiquitously observed to lie). In this case we find 
large scale magnetically-supported voids and stria- 
tions in column density imprinted on the cloud struc- 
ture (both of which have been recently observed in 
the Taurus molecular cloud), a dramatic reduction 
in the star formation rate and an effect on the IMF 
in the direction of producing fewer low mass objects. 

What effect to magnetic fields have on the dy- 
namics of the interstellar medium? Magnetic 
fields tend to smooth out ISM structure, though the 
presence of a cold gas phase can induce significant 
disorder in the global field in spiral arm shocks. 
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